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Abstract 

Recently, Gautschi introduced so-called generalized Gauss-Radau and Gauss-Lobatto for- 
mulae which are quadrature formulae of Gaussian type involving not only the values but also 
the derivatives of the function at the endpoints. In the present note we show the positivity 
of the corresponding weights; this positivity has been conjectured already by Gautschi. 

As a consequence, we establish several convergence theorems for these quadrature formu- 
lae. 
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1 Introduction 

In a recent paper [I], Gautschi considered so-called generalized Gauss-Radau and Gauss-Lobatto 
formulae which are quadrature formulae of Gaussian type, i.e., having a highest possible degree 
of exactness, and involving not only the values but also the derivatives of the function at the 
endpoints of the interval of integration. Such formulae are of the form 

/r— 1 n s— 1 

f(t)d\(t) = ^^V^W + EA^+^C-iyASi/^W + ^C/) 

3=0 3=1 3=0 

= : Qn,r,s(f) + Rn,r,s(f), (1) 
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where A is a positive measure with support being a subset of [a, b] having an infinite number of 
points of increase in (a, b), and the integers r, s > are the multiplicities of the endpoints a and 
b, respectively. In what follows we will also allow the case a = — oo (or b = +00) of a possibly 
unbounded support in which case only r = (or s = 0, respectively) is considered, that is, the 
corresponding sum in Q n ,r,s(f) does vanish. 

It is well-known and easily verified that our requirement of highest possible degree of exactness 
leads to a unique quadrature formula of the form (jTJ) with degree of exactness being equal to 
2n + r + s — 1 

V/ G P 2rl+r+s -l : Rn,r,s(f) = 0, (2) 

where here and in what follows P m denotes the space of real polynomials of degree at most m. 
Here the free nodes Tj have to be chosen as to be the simple zeros of the nth orthogonal polynomial 
with respect to the modified measure (x — a) r {b — x) s da(x) on (a, b) being clearly positive, and 
hence Tj G (a, b). Indeed, we get for Qn,o,o the classical Gaussian quadrature rule, for Q n ,i,o an d 
Qn t o,i the Gauss-Radau formulae and for Q n ,i,i the Gauss-Lobatto formula. For all these classical 
quadrature formulae it is known that the weights Ai,...,A„ and possibly Aq , , A^+j are strictly 
positive. The generalized Gauss-Radau (Gauss-Lobatto) formulae of [4] are obtained for s = 
(and r = s, respectively). 

Based on extensive numerical experiments for Jacobi, Laguerre and elliptic Chebyshev measures 
using the numerical tools and methods described in [3], Gautschi conjectured in jU Section 2.2 
and Section 3.2] that the weights of the generalized Gauss-Radau and Gauss-Lobatto formulae are 
all strictly positive. He proved himself this conjecture for the inner weights Ai, A n as well as for 
some boundary weights, namely Xq , Xq 2 ^ > for the generalized Gauss-Radau formulae Q n ,r,o, 
and Xq 1 \ Xq 2 \ A^jT 1 1 '', A„ +1 > for the generalized Gauss-Lobatto formulae Q n ,r,r- However, 
the sign of the other weights remained an open question. 

The aim of this paper is to show in Theorem [T] below that Gautschi's conjecture is true, namely, 
all weights in the quadrature formulae Q n ,r,s are strictly positive. For this we will show the slightly 
stronger result that suitable underlying Lagrange polynomials (in the Hermite sense) do not change 
sign in (a, b). As a consequence, we obtain in Corollary HJ convergence of the quadrature Q n ,r,s(f) 
for fixed r,s and n — > 00 for sufficiently differentiable functions /. The case of r,s,n — > 00 is 
discussed in Theorem [7] where we establish a geometric rate of convergence for analytic /. 

Before stating and proving our results in the next sections, we should mention that generalized 
Gauss-Radau and Gauss-Lobatto formulae are of major interest in different applications, and in 
particular in moment preserving spline approximation on a compact interval [a, b] = [0, 1], see [2], 
[31 § 3.3], and [H § 4]: given a function / G C m+1 ([0, 1]) with moments fij = J P f(t) dt, we are 
looking for a partition = r < T\ < ... < r n < r n+ \ = 1 and a spline a of class C m ~ x being 
piecewise P m on each [tj, Tj + %] for j = 0, 1, n and having the same moments 

Vj = 0,l,...,JV: / t j a(t)dt = fij, (3) 

Jo 

with N as large as possible, or in other words, the error / — a is orthogonal to Pat with respect to 
Lebesgue measure. By [3], Theorem 3.61], such a spline a exists for iV = 2n + m if and only if the 
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measure dX(t) = (—l) m+1 f( m+1 '(t)/(m\)dt on [0,1] has a generalized Gauss-Radau quadrature 
formula Q n ,m+i,m+i as in ([T]), and in this case the spline is given by the quadrature data via 

j=i j=o J - 

2 Positivity of the weights 

Theorem 1 All weights in the Gauss-type quadrature formula Q n , r , s given in (T7J) are strictly 
positive for all integers n, r, s > 



j = 1,2,.. .,n 
j = 0,l,2,...,r-l 
j = 0,l,2,...,s-l 



A, > 0, (4) 
A?> > 0, (5) 
> 0. (6) 



Proof. The property (jlj) has already been established by Gautschi [I], for the sake of complete- 
ness we repeat here the proof: consider 

= Uj {t) = (t- a y{b-ty f[ (n-tf, 

then it is clear by construction that Pjijj) = 1, and pj G F2 n +r+s-2 is non negative on (a, 6). 
According to (j2J), we may conclude that R n>rjS (pj) = 0, and thus 

Aj = Qn,r,s(Pj) = / Pj(t)d\(t) > 0, 



as claimed in (jlj). 

For a proof of (jSJ), consider the polynomial 

Pj(t) = l , ...,,.,, J U) = ^=j^nr-i-i(t)w(t), w(t) = (6 - t) s f[(r fc - t) 2 , (7) 

fc=i 

where Q m denotes the mth partial sum of the Taylor expansion of 1/u at t = a. Writing shorter 

D k f = J 

we observe that, by construction, 

Vfc=l,...,n: D^Pj = 0, Vk = 0,...,s- 1 : D%Pj = 0, VA; = 0, j - 1 : D k a Pj = 0. 
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Furthermore, for k = j,j + 1, r — 1 we find by the Leibniz product rule and by definition of f2 r 
that 



Since in addition P, G F 2n +r+s-i, we may conclude that 



A? = A? j\ DiP 3 = Q n , r , s (Pj) = | P,-(t) dA(t). (8) 
In order to discuss the sign of the expression on the right, we need the following auxiliary result. 
Lemma 2 Let P(t) = YYjLi^j with xi, ...,x m G (c, +oo) ; then 

W = 0,1,...: >0. 

Proof. For m = 1 we find that -Df (p) = ( a;i _ 1 e )i+i > 0. The general case follows by induction on 
m using the Leibniz product rule. □ 

As a consequence of the preceding lemma, we find that 

Mt)=^(t-a)^(-) 

£=0 

is strictly positive on (a, 6) for all m > 0, and thus P 3 - defined in ([7]) is also non negative in (a, b). 
It follows from (jEJ) that Aq"'' 1 > 0, as claimed in (jSJ). 

Finally, for a proof of ([6]) we observe that the variable transformation t' = —t allows to exchange 
the roles of (r, a, Aq ) and (s, b, A^) in the quadrature formula ([I]), and in particular gives a factor 
(— iy for the jth derivative. Hence the assertion ([6]) follows from (jSJ), but it is also straight forward 
to give a direct proof following the above lines. □ 

Remark 3 Notice that also D\ k Pj = for k = 1,2, ...,n. Thus, for polynomial interpolation (in 
the sense of Hermite) at the node a with multiplicity r, the nodes T\, ...,r n with multiplicity 2, 
and 6 with multiplicity s, we have shown implicitly that the Lagrange polynomials Pj = P n .r,s,a,j 
associated with the jth derivative at a do not change sign on [a, b]. It follows by symmetry that 
the Lagrange polynomial P n , r ,s,b,j associated with the jth derivative at b has constant sign (— l)- 7 
on [a, b\. However, the Lagrange polynomials associated with Tj may very well change sign on 
(a, b). 

The positivity of the quadrature weights is the essential key for proving the following conver- 
gence result both for generalized Gauss-Radau and for Gauss-Lobatto formulae. 
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Corollary 4 Let [a,b] be compact, and q := max{r — 1, s — 1}. Then for any f G C q ([a,b}) we 
have 



lim Q n , r> .(/)= / f(t)d\(t). 

n— >oo / 



Proof. In the sequal of this proof we suppose that r, s > 1, the extension of the proof for r = 
or s = is straight forward. It is not difficult to see that the space X = C q ([a, b}) equipped with 
the norm 

11/11 = max max 

0<J<<? x£[a,b] 

gives a Banach space: the completeness follows immediately from the well-known completeness of 
the space C([a, b}) with respect to the maximum norm on [a, b], see, e.g., [7J p. 258]. Also, from (TJ 
Theorem 6.3.2] it follows that polynomials are dense in (X, || • ||). Hence we are prepared to apply 
the Banach-Steinhaus Theorem: for / being a polynomial of degree k, we obtain from ([2]) that 

\/ n > k - r ~ s . g nrs(/) = J f{t) dX ( t ), 

i.e., we have convergence for a dense subset of (X, || ■ ||). For obtaining convergence in X it only 
remains to show that the norm of the linear functionals Q n ,r,s is bounded uniformly in n. Writing 
more explicitly \n), Xj{n), A^| 1 (n) for the quadrature weights occurring in (CQ), we obtain the 
simple upper bound 

r— 1 n s — 1 

\\Qn,r y s\\ < X o\ n ) + £ A i W + £ (9) 
3=0 3=1 3=0 

since, according to Theorem [H all weights occurring in these sums are positive. We observe that 

n „ 

AJV) + J2 X ^ + A «+l( n ) = ^n,r, S (l) = / dA(t) < OO. 
3=1 J 

For the remaining terms we consider the polynomial (not depending on n) 

1 — 1 s—1 

P{X) := J2 P °,r,s,a tj (x) + E i ' 1 ) ,/ 
3=1 3=1 

of degree < r+s — 1, where we recall from Remark [3] that each term in the above sums, representing 
up to a sign a Lagrange polynomial in the sense of Hermite at the abscissa a with multiplicity r 
and b with multiplicity s, is non negative on [a, b]. Hence this polynomial P is also non negative 
on [a, b], implying that, again by the positivity of the weights, 

l £^\n) + ^2x { :Un)<Q n ,rAP)= I P(t)dX(t) 

3=1 3=1 J 



for all n > 0, where in the last equality we have used §2$. Hence the expression on the right of 
is bounded uniformly in n, and the Banach-Steinhaus Theorem allows us to conclude that there 
is convergence as claimed in the assertion of Corollary HI □ 

Remark 5 By using classical arguments we may also estimate the rate of convergence of our 
quadrature formula: according to (121), we find the following bound for the error 

\Rn,rM)\ < ( [ d\{t) + \\Q n ,r,s\\) mf (10) 

V J J p&2n+r + s-l 

Here we can give a quite rough explicit upper bound for ||Q„, r)S || which is independent of n: by 
having a closer look at the construction of the polynomials Po, r ,s,a,j from (J7|) we see that 

V* e[a,b}: < P , r ,s,aAx) < E (bSS< - ~ 

and a similar bound for (— l) J Po,r,s,&,j( x )- Consequently, we learn from the previous proof and 
especially from (jHJ) that 

||Qn,r,-||< (l+r 2 (b-a) r + s 2 (b-a)^ J d\(t). (11) 

Using (TTU1) and (TTTj) . it is possible to show also convergence for a composite quadrature rule based 
on suitably shifted and scaled counterparts of Q n ,r, s , and to derive an explicit rate of convergence 
in terms of the size of the largest underlying subinterval. 



3 Rate of convergence for analytic functions 

Denote by Pf G P2 n + r +s-i the polynomial interpolating / with multiplicity r in a, multiplicity 
s in b and multiplicity 2 at the other abscissae Tj occurring in (CQ) for j = l,...,n, then using 
the Cauchy error formula for polynomial interpolation we get from ([2]) that the error for our 
quadrature formula for / G C 2n+r+s ([a, b}) may be written in terms of divided differences as 

Rn,r,s(f) = J(f(t)-Pf(t))d\(t) 

/n 
(t - a) r (t - b) s JJ(t - Tjfla^a, Tn£n, b, b, t]f d\(t) 

J =1 r 2 2s 

/n 
{t-aYit-byWit-TjfdXit) (12) 
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with £, n ,r,s,£'nrs [ a Mi since the polynomial factor in the integral is of unique sign. Denote by 
p n ,r,s the orthonormal polynomial with respect to the modified weight (t — a) r {b — t) s d\(t). We 
suppose that A has the compact support [—1, 1] with a < — 1 < 1 < b, then it is well-known from, 
e.g., 0, Section 11.11] that 

limsup|i? nA o(/)r /n <l/p 2 <l (13) 

n—*oo 

provided that / is analytic in the closed ellipse £ p with foci ±1 and half axes (p ± l/p)/2, and 
that this result is optimal for measures satisfying the Szego condition. A similar rate is shown to 
be true for fixed r,s > 0, and we are curious about the rate of convergence if r = r n and s = s n 
such that r n /n — > a > 0, s n /n — > (3 > 0. 

We first notice that £, n ,r,s £ [—1, 1] in (|T2l) . Hence, for a (set of) contour (s) C encercling once 
[—1, 1] and a, b and staying in a neighborhood of [a, b] where / is analytic, we get from the Cauchy 
formula for divided differences and from (fT2l) 



and thus 



limsup \R n ,r n ,s n (f)\ 1/n < limsupmaxr — (14) 

n^oo n^oo ^eC | Z - a \ a \ Z - bf \Pn,r n ,s n \Z) I 1 

Thus we are left with the question of nth roots asymptotics for orthogonal polynomials with 
varying weights, which has been the subject of a number of publications over the last twenty 
years, see, e.g., the monograph (8J Chapters III. 6 and VII] of Saff and Totik or the monograph [9j 
Chapter 3] of Stahl and Totik. Since the negative logarithm of the absolute value of a polynomial 
is a logarithmic potential of some discrete measure, here the right tool to describe the nth root 
asymptotic is to consider a weighted equilibrium problem in logarithmic potential theory: the 
potential and the energy of a Borel measure /j, with compact support are defined by 



U"(y)= / log( 1 :)dfi(x), J(/i) = / / log(, -)dfi(x)dfi(y). 

J F -2/1 J J F - y\ 

Define the external field Q(x) = U a (x), a = |5 a + with 5 C the Dirac unit point measure at 
x = c, then there exists a unique probability measure \x supported on S = [-1, 1] which under all 
such measures has minimal weighted energy I {pi) + 2 J Qdfi, see [U Theorem 1.1.3]. By the same 
Theorem (see also [H Theorem 1.5.1]) we also have the equilibrium conditions that U^{x) + Q{x) 
is equal to some constant F on the support of p, and > F in [—1, 1]. Then according to [5] (see 
also [9]) 

('TTS \z - aMz - W\v„ „ .AzWl) * ^{ 2 S s( U " iz) + Q(Z) ~ F ))' <15) 



with equality iff the starting measure of orthogonality A is sufficiently regular. 

For our external field, the extremal measure may be found explicitly: since Q is convex on 
[—1,1], it follows from [H Theorem IV. 1.11] that the support of p, is an interval of the form [A, B], 
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with —1<A<B<1, and also a < A and B < b since Q becomes +00 at a, b and thus in a 
neighborhood of these points we may not have equality in the equilibrium condition. 



Lemma 6 Denoting by z 1— > 0(z) = (2z - A - B + 2^(z - A)(z - B))/(B - A) the conformal 
Riemann map sending the exterior of [A, B] onto the exterior of the closed unit disk, then we have 
for z G C 



exp(2{U fl (z) + Q(z) — F)) 



<t>{z)<t>{a) 



4>{z)<j>{b) 



0(z)| 2 0(*)(0(z) - 0(a)) 0(z)(0(z) - 0(6)) 
where — 1<v4<I?<1 are defined by the system 



(16) 



a 

a 
2 1 



A - 


a 




a 


B- 


a 




a 



6- 


A 


b- 


B 


b- 


B 




A 



V. 



-V, 



= 1 ifB<l, 

<1 i/S = l, 

= 1 
< 1 



i/A>-l, 
l/^ = -l. 



Proof. Denote by a the balayage measure of a onto [A, 5], then by [SI Theorem II.4.4], fi + a 
is a positive measure of total mass (1 + f + f) supported on [A, B], and from the equilibrium 
conditions we know that its potential is constant quasi everywhere on [A, B]. However, by, e.g., [SI 
Theorem 1.1.3], the only measure satisfying this relation is (1 + | + ^)uj[a.,b], with 0J[a,b] the Robin 
measure of [A, B], i.e., the equilibrium measure with external field on [A, B\. Hence from [SI 
Eqns. (II. 4. 32) and (II. 4. 35)] and the fact that U 11 + Q — F equals zero on [A, B] we may conclude 
that 

a (3 a (3 

U^(z) +Q(z) - F = -g [AB] (z,a) + -g [A , B ](z,b) - (1 + - + -)g [AB] (z,oo), 

where by x 1— > g[A,B](x,y) we denote the Green function of the domain C \ [A, B] with pole at 
y ^ [A-S]- Taking into account the link [SI Eqn. (II. 4. 45)] between the Green function and the 
Riemann map, relation (flBT) follows. Finally, the so-called F-functional of [SI Theorem IV. 1.5] 



F(A, B) = log 



B-A, 
4 ' 



B 



Q(x) dx 



,1+ 2 + |)log( 



A \/{x 

B-A, 

— ) + 



A)(B 



a 



;9[a,b] 



a, 00 J 



g{A,B]( b i°°) 



must take its global maximum on— 1 < A < B < lat the endpoints of the support of the extremal 
measure /i. Taking partial derivatives, we arrive at the given system of equations and inequalities 
for A, B as in [SI Theorem II. 4. 4 and Lemma II. 1.15], compare with [8, Example II. 1.17] for the 
special case a = — 1 and b = 1 of Jacobi weights. □ 

In order to exploit Lemma [61 we have to consider for p > 1 the (closed) level sets £ p (a, a, 6, j3) 
being the complement of the set of z G C \ [A, B] where the right-hand side of (TIBl) is < p~ 2 < 1. 
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Figure V. Level curves for different choices of the parameters a,a } b } f3: on the left we find the 
classical case a = (3 = 0, here [A,B] = [—1,1] and the level lines are ellipses. On the right 
a = f3 = 1 and b = —a = 1, the endpoints of the support of orthogonality, here the level sets for 
p > 1 are connected with connected complement, in this case B = —A m 0.86603. 




Figure 2: Level curves for different choices of the parameters a = (3: in both cases, a = —1.5, 6 = 1, 
but on the left we find a = (3 = 0.2, on the right a = (3 = 1, leading to [A,B] = [—1, 1] in both 
cases. In particular, for small p > 1 we find three connected components for our level sets. 
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Figure 3: Level curves for different choices of the parameters (3: in both cases, a = —1.5, 6 = 1, 
a — 1, but (3 = 1 on the left (leading to [A,B] = [—1,0.8402],), and (3 = 1.2 on the right (leading 
to [A,B] = [—1, 0.79334] Here for small p > 1 we find two connected components for our level 
sets. 



Notice that for a = (3 = we have [A, B] = [—1,1], and we obtain for the complement the 
reqirement \4>(z)\ > p, that is, £ p (a, 0, 6, 0) coincides with the ellipses £ p considered before. Also, 
by the equilibrium conditions, [—1, 1] U {a, b} C S p (a, a, b, (3) for all p > 1, and from the maximum 
principle for analytic functions we may conclude that £ p (a,a,b, (3) has a connected complement 
containing a neighborhood of infinity, and at most three connected components, one of them 
containing a (if a > 0), a second 6 (if (3 > 0), and the third the interval [—1, 1], see Figure [2j 
Moreover, if A > —1 (and similarly B < 1), then p is also an extremal measure if we replace 
£ = [—1, 1] by the larger set [a, 1], and hence [a, 1] C £ p (a, a, b, (3), showing that there are only at 
most two connected components (see Figure [3]). 

If we choose as C the boundary of some level set £ p (a,a,b, (3) for some p > 1, this (set of) 
contour(s) encircles once a, b, and the interval [—1, 1]. A combination of (fT4|) . f|T5l) . and f|T6|) leads 
to the following result. 

Theorem 7 Suppose that r n /n — > a > 0, s n /n — > [3 > ; and let f be analytic in £ p (a,a,b, (3) 
for some p > 1, then 

hmsupl^^J/)! 1 /" <p~\ 

n— >oo 

One may show that again this estimate is best possible if the orthogonality measure A is 
sufficiently regular. In addition, if A satisfies the Szego condition, then following [10] we may 
obtain strong asymptotics for the orthonormal polynomials p n ,r n ,s n i an d hence with help of steepest 
descend an asymptotic equivalent of \R n ,r n ,s n {f)\- 

Different examples for the level sets of the preceding theorem are given in Figure [T], Figure [21 
and Figure [3j Though the shapes of these sets are quite different depending on the parameters, 



10 



there seem to be clearly an indication: if the function / is regular in larger neighborhoods around 
a and b, but not in such a large neighborhood around for instance (which is true for instance 
for the function f(z) = 1/(1 + z 2 )), then by the choice of larger a, (3 one improves the rate of 
geometric convergence of our generalized Gauss-Lobatto quadrature formula. 
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